Spin-wave excitations in the ferromagnetic-metallic and in the charge, orbital and 
spin ordered states in Ndi-ajSr^jMnOs with x~0.5 
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We report inelastic neutron scattering experiments on single crystals of Ndi-xSr^MnOs with 
x=0.5 and x=0.49. The spin-wave dispersion in the charge, orbital, and spin ordered state in 
Ndo.sSro.sMnOs exhibits a strongly anisotropic stiffness. The sign of the anisotropy is characteristic 
for the site-centered model for charge and orbital ordering in half-doped manganites. Within this 
model, linear spin-wave theory yields a perfect description of the experimental dispersion. In the 
ferromagnetic metallic state of Ndi-ccSra^MnOs with x=0.49 and x=0.50 magnetic excitations exhibit 
nearly the same magnon dispersion. High intense signals near the zone-boundary over a wide energy 
level overlap with a sharp spin-wave dispersion which can be described with a Heisenberg model 
including nearest-neighbor interactions. 

PACS numbers: PACS numbers: 75.10.Hk, 75.25.-j, 75.25.Dk, 75.30.Ds , 75.47.Lx, 75.40.Gb 



I. INTRODUCTION 

In spite of its eminent relevance for colossal magne- 
toresistance (CMR) the exact nature of the charge, or- 
bital and spin ordered (COSO) state in the manganites 
Ri-xAx^^Os{R=TdiTe earth, 74=bivalent alkaline earth) 
with an antiferromagnetic insulating ground state is still 
under debate. The huge change in resistivity results from 
the competition between a ferromagnetic (FM) metal- 
lic and an antiferromagnetic (AFM) insulating state. A 
comprehensive understanding of the magnetic ordering is 
essential to get a wide understanding of CMR. 

The parent compound LaMnOa with a transition- 
metal valence of Mn^+ exhibits a ground state 
which is characterized by an antiferromagnetic A-type 
ordering. ^'^ The ratio of Mn^+:Mn^+ can be controlled 
by substitution of three- valent R with two-valent A lead- 
ing to complex phase diagrams with a correlated charge, 
orbital, and spin ordered ground state (see Ref. 3 for a 
brief summary). In spite of the numerous publications 
in perovskite (113) manganites around half-doping, its 
COSO state has not been fully established.^^ Some 113- 
manganites around x=0.5 exhibit an A-type AFM metal- 
lic or FM metallic ground state while a decrease of the 
one-electron bandwidth leads to an AFM insulating state 
with CE-type ordering.^ 

The COSO state in half-doped Lai-a^Ca^^MnOa has 
been first studied by Wollan and Koehler and by Goode- 
nough proposing a simple checkerboard ordering of Mn^+ 
and Mn^+ ions (Goodenough model). ^'^^ Additionally, 
the eg orbitals on Mn^+ sites order in stripe-like patterns 
and spins order in an antiferromagnetic alignment called 
CE-type (see Fig. 1). The magnetic moments follow the 
orbital arrangement forming ferromagnetic zigzag chains 
while adjacent chains couple antiferromagnetically. First 




FIG. 1: (color online) Schematic representation of a half- 
doped perovskite manganite with charge, orbital and anti- 
ferromagnetic CE-type ordering (Goodenough model). The 
pattern shows zigzag chains propagation in [1,1,0] direction 
(cubic notation). The e^ orbital ordering on Mn^^ ions is 
shown by d33,2_^2 and d^y2_^2 orbitals. Mn^^ ions are drawn 
by balls while arrows stand for the magnetic alignment of the 
spins. The orbital unit cell is indicated by dashed black lines 
and the magnetic unit cell of Mn^^ by black dotted lines. 



structural studies confirm the Goodenough model qual- 
itatively, but the charge modulation appears to be sig- 
nificantly below the ideal value. ^'^^'^^ Single-crystal neu- 
tron diffraction experiments on Pro.6Cao.4Mn03 are not 
in agreement with the site-centered charge-ordering pic- 
ture but propose charge order on the Mn-O-Mn bonds, 
forming Zener polarons (Zener-polaron model) yield- 
ing ferromagnetically coupled dimers.^^'^^ The nature of 
the COSO state is thus still matter of debate and a 
clear proof for one of these models is still missing for 
i?i_^^^Mn03. 

The magnetic exchange coupling follows the orbital 



occupations. Consequently a careful study of these in- 
teraction parameters by analyzing the spin-wave dis- 
persion can give a direct view on the orbital arrange- 
ment. Indeed, for the single-layered (214) mangan- 
ite Lao.5Sri.5Mn04, inelastic neutron scattering exper- 
iments established the spin-wave dispersion yielding 
strong support for the COSO state described qualita- 
tively by the Goodenough model. ^^ In the magnetic 
ground state, acoustical branches exhibit a steep spin- 
wave dispersion parallel to the zigzag chains and a flatten 
dispersion perpendicular to them. These observations in- 
dicate a strong coupling inside the zigzag chains and a 
weak coupling between adjacent chains as it is expected 
for the Goodenough model. On the contrary, the Zener- 
polaron model implies a strong ferromagnetic coupling 
between the dimers and hence an anisotropy opposite to 
the observations. 

Perovskite manganites exhibit intrinsical twinning ren- 
dering microscopic investigations quite complex. The 
high-temperature structural phase transition associated 
with a lowering of the symmetry from cubic down to or- 
thorhombic results in six different domain orientations.^^ 
The occurrence of charge and orbital order can only en- 
hance the number of different domains. The magnetic 
neutron scattering at a magnetic superstructure reflec- 
tion in a twinned crystal is thus reduced by roughly a 
factor of six compared to a monodomain crystal. The 
same holds for the inelastic signal arising from a low- 
energy magnon. The interpretation of inelastic neutron 
scattering experiments is further aggravated by unavoid- 
able superposition of the signals from different domain 
orientations at higher energy. 

Due to the absence of twinning, the spin-wave 
dispersion and detailed studies in the COSO state 
has first been studied in single-layered manganites 
Lai_a^Sra;Mn04.^^'^^^^ The e^-orbitals in the undoped 
compound LaSrMn04 are dominantly aligned perpendic- 
ular to the planes, but upon doping in Lai_a^Sra^Mn04 
these e^-orbitals flop into the MnO layers^^, so that 
the ground state of the 214-manganites around half- 
doping qualitatively resembles that in 113-manganites at 
low temperatures.^^ There is thus strong support, that 
the qualitative CE-model describes the ground-state in 
single- layered half-doped Lao.5Sri.5Mn04. However it re- 
mains an open issue whether the conventional AFM CE- 
type model may account for all COSO phases in mangan- 
ites with x^O.5, or whether the order essentially changes 
with the dimensionality (three-dimensional versus two- 
dimensional materials) and with the structural distor- 
tions associated with the buckling of the MnO layers. 
Furthermore, in spite of its charge and orbital ordered 
ground-state, single-layered manganites exhibit no CMR 
under moderate magnetic field most likely due to the en- 
hanced stability of the COSO state in layered materials. ^^ 
Consequently, it is a matter of particular interest to ex- 
tend the study of the spin- wave dispersion in the COSO 
state to perovskite manganites. There is a tremendous 
number of publications on magnetic excitations in the 
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FIG. 2: (color online) Temperature-dependence of trans- 
port measurements and structural investigations of 
Ndo.sSro.sMnOs. Experimental data of electric resistiv- 
ity upon cooling show a kink at the transition into the 
(half-) metallic state at Tc=255 K (a), and the occurrence 
of the COSO state at Tco~160 K causes a sudden jump in 
resistivity by three orders of magnitude pointing to insulating 
behavior. Magnetization data indicating a ferromagnetically 
ordered state below Tc and an antiferromagnetically ordered 
state below Tco (d). X-ray powder diffraction from T=15 K 
up to T=850 K shows remarkable changes of the orthorhom- 
bic structure at various temperatures in accordance with 
electric and magnetic properties (b,c). Neutron diffraction 
at the superstructure reflection Qfm = ( 1,0,0) indicates 
ferromagnetic ordering between Tco<T<Tc. Furthermore, 
superstructure reflections referring to antiferromagnetic 
ordering of Mn^^ can be detected below Tco (e). For more 
details see text. 



ferromagnetic metallic phase and in the A-type AFM 
phase in perovskite manganites. ^'^^^^ But best to our 
knowledge, investigations of the spin- wave dispersion in 
113-manganites in the COSO state with AFM CE-type 
ordering are still missing. A microscopic analysis of both, 
the ferromagnetic metallic and the antiferromagnetic in- 
sulating state is important to understand the switching 
between theses phases by external parameters like pres- 
sure, magnetic field or temperature. 

In the paper we present experimental and theoretical 
analysis of the spin- wave dispersion in Ndi-a^Sr^^MnOa 
with x=0.50 and x=0.49. Antiferromagnetic CE-type or- 
dering in Ndi-a^Sra^MnOa has been shown to occur only 
in a very small doping range near x=0.5.^'^^'^^ A slight 
change of stoichiometry from x=0.50 to x=0.49 can lead 
to a ground-state in which ferromagnetic ordering and 
antiferromagnetic CE-type order coexist, while the fer- 
romagnetic metalhc (EMM) state above the COSO state 
is stable in a broader concentration range. ^^ The paper 
is divided as follows: After the experimental section we 
start by discussing the COSO state of the half-doped 



compound Ndo.sSro.sMnOa. We present experimental 
data of the spin- wave dispersion up to energies of about 
12 meV and compare these results with linear spin- wave 
theory. In the second part we present inelastic scatter- 
ing measurements in the FMM state of Ndo.sSro.sMnOs 
and of Ndo.5iSro.49Mn03 along the main-symmetry di- 
rections and give a description of the FM magnon dis- 
persion by a simple Heisenberg model. 



II. EXPERIMENTAL 

Two large single crystals of Ndi-a^Sra^MnOa with 
x=0.49 and x=0.50 have been grown by the floating zone 
method. The starting materials, Nd203, SrCOa, and 
Mn02 were preheated and mixed in appropriate stoichio- 
metric proportions. After sintering for several times at 
1000 °C and 1200 °C a cyhndrical rod with 90 mm length 
and 10 mm diameter was pressed and heated at 1500 °C. 
The apparatus used for crystal growth was the floating- 
zone furnace (Crystal Systems Incorporated), equipped 
with four halogen lamps. The feed and seed rods were ro- 
tated in opposite directions at about 15 rpm. The molten 
zone was vertically moved at a rate of 2-3 mm/h in a pres- 
sure of pure oxygen atmosphere (2 bar). The high quality 
of the two samples has been checked by macroscopic and 
microscopic studies. Magnetization was measured by a 
SQUID magnetometer, electric resistivity by a standard 
four-contact method, and X-ray powder diffraction ex- 
periments were carried out on a D5000 Siemens diffrac- 
tometer. 

Ndi_a,Sra,Mn03 with x=0.49 and x=0.50 crystallizes 
in the orthorhombic unit-cell with space group Ibmm 
and room-temperature lattice parameters of a^>/2xap, 
b^V2xap^ c^2xap, where ap is the lattice parame- 
ter of the cubic perovskite. In the following we use 
the cubic notations with ap=3.84 A and hence the 
reciprocal-lattice vectors are given in reduced lattice 
units as Q=(h,k,l) = |^(qa,,q^,q^). Ndi_a,Sra,Mn03 with 
x=0.49 and x=0.50 is a paramagnetic-insulator at room- 
temperature and undergoes a transition into a FMM or- 
dered state below Tc^250 K. Upon cooling a transi- 
tion of first order from the FMM state into a correlated 
charge, orbital and antiferromagnetic spin ordered state 
can be observed. ^^'^^'^^ 

All of our neutron scattering experiments have been 
performed at the Laboratoire Leon Brillouin in Saclay. 
Diffraction experiments were made using the 3T.1 diffrac- 
tometer and the G4.3 spectrometer. Inelastic neutron 
data were collected with the IT and 4F.2 spectrometers. 
The spectrometer IT, installed on the thermal channel 
of the reactor, was used with a pyrolytic graphite (PG) 
monochromator and with a PG analyzer. To suppress 
higher harmonics two PG filters were installed and the 
final neutron energy was fixed at either Ej=14.7 meV 
or E J =8.04 meV. The cold three-axis spectrometer 4F.2 
working with a PG double- monochromator was used with 
a Be-filter to suppress higher harmonics. 



III. THE CHARGE, ORBITAL AND SPIN 
ORDERED STATE IN Ndo sSro sMnOs 

The ferromagnetic metallic and antiferromagnetic in- 
sulating states in our single-crystal of Ndo.sSro.sMnOs 
have been well characterized by macroscopic techniques 
and by diffraction experiments demonstrating the high 
quality of our sample. Measurements of electric resistiv- 
ity as a function of temperature indicate a transition into 
a (half-) metallic state below Tc= 255 K and simultane- 
ously a rapid increase of magnetization points to ferro- 
magnetic ordering [Fig. 2(a),(d)]. The transition into 
a charge, orbital and spin ordered state at Tco^lSO K 
is of the first order. Herby, the electrical resistivity ex- 
hibits a jump by more than three orders of magnetiude 
and the magnetization vanishes. X-ray powder diffrac- 
tion has been performed from T=15 K up to T=850 K. 
The refinement of the data was carried out with the Ri- 
etveld method implemented in the FullProf suite. ^^ The 
best description of the data can be obtained with the 
space group Ibmm for the studied temperature range. 
The symmetry reduction from the space group Pm3m 
of an ideal undistorted perovskite to space group Ibmm 
arises from tilting of MnOe octahedra around an axis 
parallel to the edges in the ab-planes. The orthorhom- 
bic distortion of the lattice parameters increases rapidly 
at T=550 K which can arise from an orbital ordering 
effect of Mn^+ ions [Fig. 2(c)]. Upon cooling through 
the FMM phase transition, a slight increase of the or- 
thorhombic c-axis with a shrinking of the a and b-axis is 
observed. It appears likely that upon entering the FMM 
state an enhancement of the in-plane metallicity an FM 
exchange is favored. An abrupt change of the lattice pa- 
rameters during the transition into the COSO state indi- 
cates a sudden and strong structural change [Fig. 2(b)]. 
Temperature-dependent neutron diffraction experiments 
at the Bragg reflection Qfm=( 1,0,0) prove ferromagnetic 
ordering between Tco<T<Tc which is in perfect agree- 
ment with our magnetization data [Fig. 2(e)]. Note, 
that the FM signal at Qfm = ( 1,0,0) is also visible above 
Tc, but it gets strongly suppressed in the COSO state. 
The temperature-dependent behavior of the superstruc- 
ture reflection intensity of Q7Vfn3+ = (0.25,0.25,0.5) pro- 
vides information on magnetic ordering of Mn^+ below 
Too [Fig- 2(e); for more details see Ref. 19]. The quality 
of our sample is evident from the observed sharp transi- 
tions, the increase in resistivity by several orders of mag- 
nitude, and from the perfect agreements of Tc and Too 
with pubHshed data.^^'^^ Although charge and orbital or- 
dering cannot be directly detected by neutron diffraction, 
charge and orbital ordering leads to a characteristic dis- 
tortions of the MnOe-octahedra that are visible in neu- 
tron diffraction experiments due to the higher sensitivity 
to oxygen. 

The Goodenough model predicts an alternating charge 
ordering of Mn^+ and Mn^+ ions. Considering an ini- 
tially ideal cubic perovskite, this doubles the structural 
unit cell to {V2ap x ^/2ap x ap) and causes superstruc- 
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FIG. 3: (color online) Dispersion of the magnetic excitations in Ndo.sSro.sMnOs along the main- symmetry directions (a). The 
path F-A stands for investigations in [1,0,0] direction (diagonal to the chains), while F-B measures the dispersion perpendicular 
and F-C parallel to the zigzag chains. Investigations parallel to the c-axis are presented by F-Z. Solid curves are fits as 
described in the text and crystal field excitations of the Nd ions are indicated by grey-shaded bars. The different paths along 
the main-symmetry directions around an antiferromagnetic zone-center for chains propagation in [1,1,0] directions is shown in 

(b). 



ture reflections with propagation vector kco=(±^,±^,0). 
Crystallographic studies find far smaller disproportion 
than a complete Mn^+/Mn'^+ (see Ref. 11,12). Neverthe- 
less we use the nominal labels Mn^+ and Mn^+ through- 
out this paper. The orbital ordering of e^-electrons on 
Mn^+-sites reduces the symmetry to orthorhombic with 
lattice parameters 2>/2ap along [1,1,0] and >/2ap along 
[-1,1,0] for zigzag chains running along [1,1,0]. As the 
zigzag chains can propagate along [1,1,0] and [-1,1,0], 
orbital superstructure reflections with propagation vec- 
tor koo = (i|,i|,0) can be detected. The orthorhombic 
distortion induces a twofold twinning since the orbitals 
can stack either along [1,1,0] or [-1,1,0] direction. Both 
arrangements superimpose and contribute to an equal 
amount in the scattering geometry studied in our sam- 
ple. The superposition with the structural twinning in 
these pseudo-cubic compounds can lead to a very com- 
plex twinned structure in the COSO state. The charge 
and orbital ordering implies an antiferromagnetic CE- 
type ordering with ferromagnetic alignment of magnetic 
moments along the chains and antiferromagnetic cou- 
pling of adjacent chains. Fig. 1 shows a schematic 
representation of the charge, orbital and spin ordering 
with zigzag chains propagation in [1,1,0] direction. In 
perovskite manganites the planes of the orbital order- 
ing always align with the structural distortions due to 
the octahedron tilting. The latter leads to an unit cell 
V^ttp X \f2a^ X 2ap with the doubling being always per- 
pendicular to the layers of orbital order (note that the 
zigzag chains run in these). Due to the strong struc- 
tural distortion, orbitals and charges align in rows along 
this perpendicular direction (c-axis). Charge and orbital 
order are thus not modulated along this direction, cor- 
responding to a ferro-type stacking. In the following the 
stacking direction is taken always along [0,0,1] in our cu- 
bic notation. ^^ The magnetic CE-type spin ordering re- 



sults in two sublattices: the magnetic unit cell of Mn^+ 
and the magnetic sublattice of Mn^+. The latter one is 
depicted in Fig. 1. Superstructure reflections referring to 
the spin ordering of Mn^+ can be detected at propagation 
vectors kMn4+=±(^,0,^) or \(.y^^^^=^(^,\,^ and super- 
structure reflections referring to magnetic ordering of 
Mn^+ at kMn3+=i(|7i|7^)- To investigate pure charge 
or orbital ordering we have to measure superstructure re- 
flections with large Q where magnetic order is negligible 
because the magnetic form-factor decreases rapidly with 
increasing scattering vector. 

To study the spin- wave dispersion in the AFM ordered 
state along the main-symmetry directions, we aligned 
the sample in two different scattering-planes. Most of 
our experiments have been done by aligning the crystal 
in the (1,1,0)/(0,0,1) scattering plane to access momen- 
tum transfers of the kind (h,h,l). In order to focus on 
the spin- wave dispersion parallel to the zigzag chains the 
sample was orientated in a less conventional geometry 
with reflections (2, 1,1)/ (0,1,-1) defining the scattering 
plane. The spin- wave dispersion at T=12 K, observed 
by inelastic neutron scattering from 2 meV up to about 
12 meV is depicted in Fig. 3. The solid lines corre- 
spond to a fit based on Hnear spin-wave theory. Details 
on the theoretical calculations of the magnetic excitation 
spectrum are given below. Besides the superposition of 
all twins rendering the study at high energies impossi- 
ble, we observe Q-independent signals around 0.9 meV 
and around 12 meV pointing to crystalline electric field 
(CEF) excitations of the Nd ions. Scans in different 
Brillouin zones and observations in similar compounds 
support this conclusion. ^^'^^'^^ Consequently clear mag- 
netic excitations in all symmetry directions at 12 meV 
are not recordable. Our spin- wave studies have been 
performed around several antiferromagnetic zone centers, 
F = (|, |, ±^) with h,k odd integer, and have been com- 
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FIG. 4: (color online) Representative energy and constant- 
energy scans at T=12 K in Ndo.sSro.sMnOs along the main- 
symmetry directions. Two constant- Q scans which measures 
the dispersion along F-B are depicted in (a). Constant-energy 
scans to study the acoustic branch parallel to the chains (F- 
C), 45° diagonal to the chains (F-A) and in c-direction (F- 
Z) are shown in (b)-(c). Solid lines correspond to fits with 
Gaussians. 



bined to the presented dispersion. Fig. 4 shows examples 
of typical raw-data scans to study the spin- wave disper- 
sion along the main- symmetry directions. We were able 
to separate between different paths since only one of the 
orbital twins contributes to a quarter-indexed magnetic 
zone center. For example, to measure the spin- wave per- 
pendicular to the zigzag chains (F-B) starting at the zone 
center (^, |, ^) one has to measure along the [-1,1,0] di- 
rection while the [1,1,0] direction measures the path par- 
allel to the zigzag chains (F-C) [Fig. 3(b)]. Starting from 
the antiferromagnetic zone center, we find a gap and a 
splitting into two single modes, similar to observations 
in the single-layered system Lao.5Sri.5Mn04 (Ref. 15). 
The degeneracy of the two AFM magnon branches seems 
to be removed due to the single-ion anisotropy. Energy 
scans at several magnetic zone-centers exhibit two well 
distinguishable contributions at energies Ei=2.4(l) meV 
and E2=4.1(l) meV [Fig. 5(a)]. The gap and energy 
splitting of the lowest excitations can be explained by the 
alignment of the magnetic moments. The easy-axis of the 
magnetic moments of Ndo.sSro.sMnOa is in the ab-plane 
at an angle of 45° with respect to the a and b axis (cubic 
notations).^^ To study the splitting near the zone-center 
and to follow the two modes in the directions F-B and 
F-C we performed inelastic neutron experiments with the 
cold-source spectrometer. Away from the zone-center, at 
q>0.025, the two excitations merge into a single signal. 

The lowest acoustic branch along the path F-C exhibits 
a steeper dispersion than that along the path F-B. The 
endpoint of the acoustic branch at B can be observed at 
7.5(1) meV. At this energy, we find the acoustic branch 
along F-C in the middle of the Brillouin zone pointing to 



a stronger ferromagnetic coupling parallel to the zigzag 
chains in comparison to the coupling perpendicular to 
the chains. Although magnetic excitations at the zone- 
boundary parallel to the chains can experimentally not be 
resolved due to the twinning, the anisotropic behavior is 
apparent. Our theoretical calculations predict a magnon 
energy of about 22 meV at the zone boundary which is 
three times larger than the endpoint B of the branch F-B. 
This study is in full agreement with the observations in 
Lao.5Sri.5Mn04 (Ref. 15). The antiferromagnetic cou- 
pling along the c-direction (F-Z) exhibits an intermedi- 
ate dispersion: this means a steeper acoustic branch in 
comparison to the antiferromagnetic coupling in F-B, but 
flatten than the ferromagnetic coupling along the zigzag 
chains (F-C). Consequently the antiferromagnetic cou- 
pHng between adjacent chains is weaker in comparison to 
the antiferromagnetic stacking in c-direction. This is fun- 
damentally different in the layered materials, where there 
is no significant c-axis coupling. The CE-type order in 
Ndo.sSro.sMnOs is thus - as expected - much more three- 
dimensional in nature, the AFM stacking in c-direction 
is stronger than the AFM coupling between zigzag chains 
in the ab-layers. 

Theoretically, the formation of the CE-type super- 
structure close to half filling is understood in terms of 
the degenerate double-exchange (DDEX) model. ^'^^'^^'^^ 
The t2g electrons form local magnetic moments (core 
spins) which are coupled to the itinerant electrons of 
the degenerate eg orbitals by a ferromagnetic interatomic 
Hund's exchange Jh- In the context of the mangan- 
ites, the DDEX model can be greatly simplified since 
Jh is much larger than the electronic bandwidth and the 
super-exchange between core spins which are treated as 
classical. In the limit Jh -^ oo, the itinerant-electron 
spins are perfectly aligned with the local- moment spins. 
Therefore, the itinerant electrons can move only in a fer- 
romagnetic background. The complex CE superstructure 
is the result of the frustration between the electronic ki- 
netic energy and the antiferromagnetic super-exchange 
between the core spins. In this phase, the electronic ki- 
netic energy is directed along ferromagnetic zigzag chains 
which form because of a cooperative ordering of charge- 
and orbital degrees of freedom. ^^ 

Although, the DDEX model treated in the semiclassi- 
cal approximation can explain the richness of the phase 
diagram^^ and in particular account for the CE struc- 
ture observed in the close vicinity of half filling, ^'^^'^^'^^ 
the calculation of the magnetic excitation spectra within 
the DDEX model is extremely difficult since it requires 
finite Jh and quantum mechanical treatment of the core 
spins. ^^'^^ To avoid these comphcations, we compute 
the spectrum from a spin-only model with effective ex- 
change couplings which includes both, AFM-super- and 
FM-double-exchange contributions. The latter are ex- 
pected to contribute predominantly along the chains and 
to over compensate the bare AFM super-exchange. Con- 
sequently, a minimal effective model consists of FM cou- 
pHngs Jfm along the chains and AFM coupHngs Jafm 
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FIG. 5: (color online) Energy scans in the first Brillouin zone 
at the antiferromagnetic center Q= (0.75,0. 75,0. 5) (a) and at 
the zone boundary B along the acoustic branch F-B (b). The 
gap and the splitting of the excitations at the zone-center is 
clearly observable (indicated by black arrows). The signal 
around E=l meV corresponds to CEF excitations of the Nd 
ions. (The elastic line is not shown). The excitation at the 
end-point B can be observed at 7.5(1) meV while the signal 
around 12 meV can be ascribed to a CEF, as well (b). Dashed 
lines correspond to individual fits with Gaussians and the 
solid line to the resulting fit. 



between the chains. 

Indeed, the magnetic excitation spectrum of the lay- 
ered system Lao.5Sri.5Mn04 can be explained remark- 
ably well by this minimal model. ^^ Including an ad- 
ditional next-nearest neighbor exchange Jfm,2 between 
Mn^+ sites along the chains (Fig. 6) and a small single- 
ion anisotropy A to account for the spin gap, one finds an 
excellent description of the neutron scattering data even 
for optical branches at higher energies which has been 
interpreted as further evidence for the existence of the 
AFM CE-type superstructure at half filling. 

In the following, we use exactly the same model to fit 
the spectrum of the perovskite Ndo.sSro.sMnOs including 
an additional coupling Jafm,c along the c direction (Fig. 
6). The effective spin-only Hamiltonian reads 



+ Jafm 2_^ ^^^J + Jafm,c 2_^ ^*^i 

+ AE(^f)' (1) 

i 

where the exchange couplings Jfm^ Jfm, 2^ Jafm^ 
Jafm,c and single- ion anisotropy A are defined to be 
positive. The sites on the three dimensional lattice are 
labelled by z, (•, •)|| and ((•, •))|| denote nearest and next- 
nearest neighbor bonds along the FM zigzag chains. Note 
that Jfm,2 acts only between Mn^+ spins on the corners 
of the chains. Nearest-neighbor bonds between the chains 
in the a6-plane and along the c direction are denoted by 
(•, •)^ and (•, •)c, respectively. In order to facilitate the 
comparison with Ref. 15 we use the same definition of 
the sum: each bond appears only ones. 

In the following, we use linear spin- wave theory to 



calculate the magnetic excitation spectrum and subse- 
quently fit the experimental data. The procedure is stan- 
dard and details can be found elsewhere. ^^'^^ After repre- 
senting the spins by Holstein-Primakoff bosons the spec- 
trum is obtained by a Bogoliubov transformation which 
amounts to diagonalizing a 16 by 16 matrix which is ex- 
plicitly given in Ref. [45] in the case Jfm, 2 = A = 0. 

Whereas the idealized CE structure assumes a checker- 
board ordering of Mn^+ and Mn^+ sites with correspond- 
ing local moment spins of 6'i = 2 and S2 = 3/2, respec- 
tively, the experimentally observed charge and spin mod- 
ulations are significantly smaller. To account for this, in 
the following we use Si = 1.79 and 5*2 = 1.67 such that 
5*1 5*2 = 3 and S1/S2 = 3/2.8 equal to the ratio of the 
measured moments. ^^'^^ 

The best fit parameters are summarized in Table I. 
As in the layered system, ^^ the FM coupling along the 
chains Jfm is found to be much larger than the AFM 
in-plane coupling Jafm between the chains. For the 
Lao.5Sri.5Mn04 compound we used for the spins of Mn^+ 
Si = 2 and for the spins of Mn'^+ ^'2 = 1.5. As the 
product of 5'i5'2 = 3 is the same for both systems, the 
nearest-neighbor parameters Jfm and Jafm can directly 
be compared. 

This observation and the ratio Jafm,c/Jafm ^ 4 
are in agreement with a quantitative estimate of the ex- 
change couplings based on a modified Kanomori model 
which carefully takes details of the lattice structures of 
various CE-type manganites into account. ^^ 





Ndo. 


sSro.sMnOs 


Lao. 


5Sri.5Mn04 


Jfm [meV] 




6.76 




9.98 


Jfm,2 [meV] 




2.76 




1.83 


Jafm [meV] 




0.75 




3.69 


Jafm,c [meV] 




3.22 




- 



TABLE I: Results of the obtained fit parameters in the AFM 
CE-type state for Ndo.sSro.sMnOs in comparison to the fit 
parameters in Lao.5Sri.5Mn04 from Ref. 15. 

Jafm,c connects either two Mn^+ or two Mn^+ sites, 
whereas Jafm couples one Mn^+ and one Mn^+ site 
where imperfect orbital polarization will imply a FM 
double-exchange interaction rendering the AFM super- 
exchange. As in addition the crystal lattice is consid- 
erably flattened, Jafm,c is much larger than Jafm- 
In our model we include a simple single-ion anisotropy 
term A^^(S'f)^. The ordered moments in the CE-type 
phase in Ndo.sSro.sMnOs are aligned along the ab-plane. 
Therefore an anisotropy term perfectly accounts for an 
ab-easy-plane system. However the spins seem to be well 
fixed also within the ab-layers causing a splitting in the 
zone-center magnon energies. ^^ The same situation is also 
found in single-layered Lao.5Sri.5Mn04, where the lower 
excitation could be assigned to the in-plane anisotropy. 
Also in Ndo.sSro.sMnOs we may identity the lower exci- 
tation with the in-plane anisotropy. It costs more energy 
to rotate the spins into the c-direction. A full treatment 
for such complex anisotropy has been given for a simple 
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FIG. 6: (color online) Illustration of the effective spin-only 
model for the AFM CE-type structure. Spins along the zigzag 
chains are coupled by ferromagnetic nearest and next-nearest 
neighbor exchanges Jfm and Jfm,2, the latter acting only 
between Mn^^ sites. The chains are coupled antiferromag- 
netically by exchanges Jafm and Jafm,c between neighbors 
in the a6-plane and along the c-axis, respectively. 



AFM layered structure in Ref. 48. 

Inelastic neutron scattering experiment have 
been able to determine the spin-wave excitations of 
Ndo.sSro.sMnOs for energies up to about 12 meV. Our 
data exhibit a steeper magnetic dispersion parallel to the 
zigzag chains while the dispersion perpendicular to the 
chains is flatten. This observation is in contrast to the 
Zener-polaron model which implies a strong coupling in- 
side the dimers. Taking this and the magnetic structure 
of the Zener-polaron scenario into account one obtains 
just the opposite ansisotropy in the dispersion slopes 
along the F-B and F-C directions studied here. The 
spin-wave theory within the CE-type structure explains 
perfectly our experimental data. The experimental 
and the theoretical analysis are in agreement with the 
single-layered compound Lao.5Sri.5Mn04 (Ref. 15). 
Consequently we may extend the strong support of the 
Goodenough model to perovskite manganites around 
half-doping. 



IV. THE FERROMAGNETIC METALLIC 
STATE IN Ndo.5oSro.5oMn03 

In order to analyze the competition between the insu- 
lating COSO and the FMM phases in Ndo.sSro.sMnOs 
we performed studies of the magnon dispersion in the 
FMM phase as well. The magnetic excitations in the 
FMM state in perovskite manganites have been stud- 
ied in many different compounds and a damping of the 
spin-wave dispersion towards the zone boundaries have 
been reported in many cases. ^^'^^'^^^^'^^'^^ To explain 
this behavior, several microscopic mechanisms have been 
proposed, without reaching a clear conclusion (for a brief 
summary see Ref. 51). One explanation for the soften- 
ing of the spin- wave dispersion might be the development 
of quantized spin- waves within ferromagnetic clusters of 
restricted sizes as reported in Ref. 27,31. 
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FIG. 7: (color online) Full spin- wave dispersion along 
the main-symmetry directions in the FMM state in 
Ndo.sSro.sMnOs at T=174 K. Experimental data were taken 
along the directions F-R, F-M, and F-X, respectively. Solid 
lines refer to a fit with a Heisenberg ferromagnet with only 
nearest-neighbor interactions and dotted lines to a fit with 
nearest and forth-nearest-neighbor interactions. A significant 
improvement of the fit to the data by taking additionally J4 
into account can not be observed. An enhanced signal around 
the zone boundary might be given by several magnetic levels 
(b). For more details see text. 



We studied the spin dynamics of Ndo.sSro.sMnOs at 
T=174 K, where the ferromagnetic signal of the reflection 
at Qfm=(1,0,0) exhibits its maximum [see Fig. 2(e)]. 
To investigate magnetic excitations along the main- 
symmetry directions [1,0,0], [0,1,1] and [1,1,1], we aligned 
the sample with reflections (1,0,0) and (0,1,1) defining 
the scattering plane. All data have been collected around 
the Bragg-reflections Qfm=(1,0,0), Qi?M=(0,l,l) and 
Qfm= (1,1,1)- Most data have been collected with the 
thermal spectrometer IT with fixed Ej=14.7 meV. The 
calibration of the spectrometer has been checked by 
a Vanadium standard crystal. An energy scan with 
Ej=14.7 meV indicates a resolution of AE^l meV. Fig. 
7 presents the spin- wave dispersion along the main sym- 
metry directions. The data have been corrected by the 
Bose factor and finally fitted with the ResLib code which 
calculates the convolution of the four-dimensional reso- 
lution function with the dispersion. ^^ We were able to 
investigate magnetic excitations of the acoustic branch 
along the pathes F-R, F-M, and F-X. An energy scan 
with the cold-source spectrometer 4F.2 at the zone-center 
Qfm=(1,0,0) yields no evidence for an excitation gap. 
Constant-energy scans around the ferromagnetic zone- 
center Qfm = (1,0,0) are depicted in Fig. 8(a). Investi- 
gations in different Brillouin zones clearly identify a mag- 
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FIG. 8: (color online) Constant-energy scans around the zone- 
center Qfm = (1,0,0) in Ndo.sSro.sMnOs along the symmetry 
direction [1,0,0]. Data at E=l meV and 3 meV are recorded 
with the spectrometer IT with fixed E/=8.05 meV and at 
E=8 meV with fixed E/ = 14.7 meV, respectively (a). A de- 
crease of the signal with increasing scattering vector demon- 
strates its magnetically origin (b). 



h = -J2j^ 



I J O^Oj 



(2) 



Herby Jij describes the magnetic exchange coupling be- 
tween two spins Si and Sj at the position R^ and Rj. 
Note that a pair appears twice in the sum.^^ In linear 
approximation the spin- wave dispersion with only ferro- 
magnetic coupling {Jij > 0) becomes: 



^(q) = nu{q) = A + 2S[J{0) - J(q)] 



(3) 



with: 



-/(q) = T.-^^^ 



iq(Ri— Rj 



netic origin for the signals appearing at low energy trans- 
fer, since the intensity follows the magnetic form factor 
as shown in Fig. 8(b). In Ref. 25 the spin- wave disper- 
sion in the FMM state in Ndo.sSro.sMnOs along different 
symmetry directions has been studied and an extremely 
flat dispersion with energies around 10 meV at the zone 
boundaries is reported. Our experiments disagree with 
theses observations since our data yield much steeper 
branches. For example, clear excitations along the [1,1,1] 
direction between energies of E=18 meV and E=22 meV 
can be observed in the middle of the Brillouin zone (see 
Fig. 9). The main magnetic signal seems to stem from 
the simple ferromagnetic magnon branch; but there are 
weak but significant additional signals which are most 
likely magnetic as well. The constant energy cuts (h,0,0) 
across the (1,0,0) FM Bragg-peak show more intensity 
near the magnetic zone boundaries ^(0.5,0,0) than at 
the center which can not be attributed to the back- 
ground. Also the energy scan at (1.5,0,0) exhibits sig- 
nificant intensity below the energy of the magnon branch 
[Fig. 7 (b)]. This scattering resembles the splitting of 
the magnon branches observed in Lai_a;(Ca,Sr)a;Mn03 at 
lower doping, which was attributed to a confinement. ^^'^^ 
However, such confinement effect - if present - must be 
much weaker in this higher-doped Ndo.sSro.sMnOs ma- 
terial. 

Spin- wave dispersion in 113-manganites with a large 
electron bandwidth and high Tc can be described by cal- 
culations with the conventional Heisenberg model includ- 
ing only the nearest-neighbor exchange inter actions. ^^'^^ 
In the following the description of our data takes only 
the sharp signals into account without keeping the in- 
creasing intensity near the zone-boundaries in mind. To 
evaluate the spin- wave dispersion in Ndo.sSro.sMnOa we 
use two models: (1) we consider only the ferromag- 
netic coupling between nearest-neighbor interaction Ji, 
(2) nearest-neighbor interactions Ji and forth-nearest in- 
teractions J4 are taken into account. 

The Hamiltonian of a Heisenberg ferromagnet can be 
written as: 



A finite energy gap near the zone-center can be con- 
trolled by A. Since we do not observe a finite gap for the 
[1,0,0], [0,1,1] and [1,1,1] directions, we will set A equal 
to zero. 

The simplest case with only nearest-neighbor Ji inter- 
actions leads to the relation: 



E{q) = A^4SJi[3-cos{qha)-cos{qka)-cos{qia)] (4) 

Taking only Ji into account, we can very well describe 
our experimental data (Fig. 7 solid lines). The fits of the 
data with only the nearest-neighbor interaction Ji yield 
26'Ji=7.6(2) meV for the [1,0,0], 2SJi=6.3{2) meV for 
the [0,1,1] and 26'Ji=7.2(2) meV for the [1,1,1] direction. 
To compare an improvement of the description of the fit 
by introducing the coupling J4 as suggested previously 
in Ref. 50, the dispersion for the [1,0,0] direction takes 
the form [Fig. 11(b)]: 



^(q) 



45Ji[(l 

45J4[(1 



cos{27rqh))] + 
cos{4.7rqh))] 



(5) 



The best fit to the data with nearest-neighbor and 
forth-nearest interaction interaction is depicted in Fig. 
7 by the dotted lines. There is only little difference in 
the spin- wave dispersion of the two models in the energy 
range studied. Therefore our data give no support for a 
related 2S'J4 parameter. 

In the CE-type state we obtained the best fit to the 
data for the ferromagnetic coupling parallel to the chains 
with 26'Jfm=11-7 meV. In the FM state we can describe 
the data with a mean exchange energy of 2571=7.0(5) 
meV for the fit with nearest-neighbor interaction tak- 
ing into account. This means the nearest-neighbor ferro- 
magnetic exchange in the FM state is strongly reduced 
by a factor of about 1.5 in comparison to the AFM 
CE-type state, but note that each spin has six ferro- 
magnetically coupled spins in the FMM state, compared 



1000 



750 - 






500 - 



c 

CD 

c 250 




0.05 0.1 0.15 0.2 0.25 
h [r.l.u.] 

FIG. 9: (color online) Observed excitation at high energies 
along the symmetry direction [1,1,1]. The maximum shifts 
to increasing q with increasing energy indicating the steeper 
dispersion. 



to only two FM neighbors in the CE-phases. In the 
single-layered compounds Lao.5Sri.5Mn04 an exchange 
energy of 2S'Ji?M=7.5(5) meV has been observed above 
the AFM CE-type ordered state. ^^ Consequently, the ex- 
change parameters in the FM states in Lao.5Sri.5Mn04 
and in Ndo.sSro.sMnOa are almost identical. This re- 
sult is quite interesting since the single-layered compound 
does not show metallic behavior and possesses less neigh- 
bors. More importantly, the layered compound only ex- 
hibits ferromagnetic fluctuations generated by ferromag- 
netic clusters with finite size above Tco- Our inves- 
tigation in the FMM state in Ndo.sSro.sMnOs indicate 
a quite good description of the data along the main- 
symmetry directions by taking only the nearest-neighbor 
coupling Ji into account (but note that our model is 
slightly over-simplified since we do not take the inten- 
sities near the zone-boundary taking into account). 

To get a more precise perception of the magnetic exci- 
tation in the FMM state in Ndi-a^Sra^MnOa and to com- 
pare the exchange parameters we decided to study the 
spin-wave dispersion in a slightly electron-doped single 
crystal with a composition of Ndo.5iSro.49Mn03. 



V. THE FERROMAGNETIC METALLIC STATE 

IN Ndo.5iSro.49Mn03 

A single-crystal of Ndi-a^Sr^^MnOs with doping 
level x=0.49 has been characterized by temperature- 
dependent measurements of magnetization, resistivity. 
X-ray, neutron diffraction and by EDX analysis. A sig- 
nificant drop of resistivity accompanied by an increase of 
magnetization reminiscent of the FMM state is observed 
at Tc?^245 K (Fig. 10). In accordance with the macro- 
scopic studies, the scattering at Qi?M=( 1,1,0) exhibits 
an increase of intensity below Tc which can be ascribed 
to ferromagnetic correlations [Fig. 10(c)]. Upon further 
cooling, the electric resistivity increases by more than 
two orders of magnitude at the first order transition into 
the COSO state at Tco^ll2 K. Temperature-dependent 
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FIG. 10: (color online) Characterization of the single crystal 
Ndo.5iSro.49Mn03 by macroscopic and microscopic studies. 
A ferromagnetic metallic state can be observed by electric re- 
sistivity measurements and magnetization data (a,c). Elastic 
neutron scattering at Qfm = ( 1,1,0) confirms ferromagnetic 
ordering below Tc (d). A sudden increase in intensity of 
the superstructure reflection Qoo= (2.25,0.25,0) below Tco 
refers to orbital ordering (b),(d). 
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FIG. 11: (color online) Dispersion of the magnon excitation 
in the FMM state in Ndo.5iSro.49Mn03 along the symme- 
try directions [1,0,0] and [0,1,1] at T=151 K (a). Solid lines 
refer to fits with only nearest-neighbor interactions Ji and 
dashed lines point to calculations with nearest and forth- 
nearest-neighbor interactions Ji and J4 as depicted in (b). 
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FIG. 12: Raw-data scans with fixed energy observed in [1,0,0] 
direction. Data at 2 meV are taken with E/=8.05 meV and 
above with E/ = 14.7 meV. Intensities are normahzed to same 
monitor but shifted as guide to the eyes (a). Signals along 
the symmetry direction [0,1,1] exhibit a loss in intensity and 
increase of width with increasing energy, too (b) . 



investigations of the superstructure reflection referring to 
orbital ordering Qoo = (2. 25, 0.25,0) yield a sizeable sig- 
nal below Too [Fig. 10(d)]. Comparing the transition 
temperatures into the COSO state in Ndi-a^Sra^MnOa for 
x=0.50 and x=0.49, we observe a less robust COSO state 
in Ndo.5iSro.49Mn03. Before discussing inelastic data we 
want to compare the ground state properties of the ide- 
ally doped system (x=0.5) in comparison to the slightly 
electron-doped system (x=0.49). In Ndo.5iSro.49Mn03 
below T!^50 K, the electric resistivity is almost constant. 
In contrast, for x=0.5 we find a continuous increase of re- 
sistivity below Tco up to p~10^ Qm at the lowest tem- 
perature in comparison to p^lO^ ftm. for x=0.49. The 
slight doping difference between these two compounds 
strongly affects the magnetization at low temperatures 
as well. The magnetization with an applied field of B=50 
mT saturates with m=0.05 /i^/Mn for x=0.49 while the 
magnetization for x=0.50 tends to zero at low tempera- 
tures. These properties underline the coexistence of anti- 
ferromagnetic CE-type and ferromagnetic ordering in the 
ground state of Ndo.5iSro.49Mn03 as reported in Ref. 32. 
Neutron diffraction of the reflection Qi?M=( 1,1,0) fur- 
ther supports this interpretation. The signal below Tc 
increases as in the half-doped system, but in the COSO 
state in x=0.49 a high intensity can still be recorded 
at low temperature. We analyzed the stoichiometry of 
the samples with EDX analysis. However the accuracy 
of this method does not allow for determining the con- 
tent of oxygen concentration. Hence we may not rule 
out an oxygen- deficiency. Concerning the various tran- 
sition temperatures and the temperature dependence of 
the electrical resistivity our sample is fully comparable 
to the Ndo.5iSro.49Mn03 crystal studied in Ref. 32. 

We have studied the spin-wave excitations in 
Ndo.5iSro.49Mn03 at T=151 K, which is well in the ferro- 
magnetically ordered state, and present an analysis of the 
data with the same Heisenberg models as those used for 
the half-doped compound. The scattering experiments 
have been performed by aligning the crystal with the re- 
flections (1,0,0)7(0,1,1), allowing momentum transfers of 



the kind (h,l,l). We focused our measurement on exci- 
tations along the symmetry directions [1,0,0] and [0,1,1] 
in cubic notation. Data have been measured in the Bril- 
louin zones around the zone-centers Qfm = (1,0,0) and 
Qfm= (0,1,1) and were finally combined. The data have 
been corrected by the Bose factor and finally fitted with 
the convolution of the four-dimensional resolution func- 
tion with the dispersion. Fits to the data results in the 
spin- wave dispersion as shown in Fig. 11. Raw-data 
scans with fixed energy around ferromagnetic zone cen- 
ters along the [1,0,0] and [0,1,1] directions are depicted 
in Fig. 12. 

In the following we will proof the experimental data 
of the sharp signals with the Heisenberg model as de- 
scribed for the half-doped system. The investigations 
of the dispersion near the zone center F with the cold 
triple-axis spectrometer are in agreement with our ob- 
servations in Ndo.5Sro.5Mn03 and we find no indication 
for a spin- wave gap in Ndo.5iSro.49Mn03. Consequently, 
we set A equal to zero in our calculation. To evaluate 
the spin- wave dispersion in Ndo.5iSro.49Mn03 we make 
the same assumptions as those used for the x=0.5 com- 
pound: (1) we consider only the ferromagnetic coupling 
between nearest-neighbor interaction Ji, and (2) nearest- 
neighbor interactions Ji and forth-nearest interactions J4 
are taken into account. 

All fit parameters for the spin-wave dispersion in 
Ndo.5iSro.49Mn03 and Ndo.5Sro.5Mn03 are summarized 
in Table II. Since no significant improvement of the re- 
duced chi-squared Xred ^^^ ^^^ model including J4 and 
since the ratio of J4/J1 is quite small, there is no evidence 
that the forth-nearest-neighbors is needed to describe the 
ferromagnetic spin- wave dispersion in Ndi-^^Sr^^MnOs 
with x=0.49. 
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2SJi [meV] 
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7.6(2) 
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model 2 

2SJi [meV] 

2SJ4 [meV] 
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6.0(3) 

0.5(2) 

8 
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0.8(1) 

16 


8.0(3) 

-0.2(2) 
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5.3(9) 
0.4(4) 
8 


2.1(9) 
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19 



TABLE II: Fit parameters for the spin- wave dispersion in 
Ndi_a,Sra;Mn03 with x=0.49 and x=0.50 in the EMM state. 
The data have been obtained from two different models. 
Model 1, only nearest-neighbor interactions Ji and model 
2, nearest-neighbor interactions Ji and forth- nearest inter- 
actions J4 are taken into account. 

Our studies for x=0.49 yield the identical dispersion 
in comparison to the optimal-doped compound x=0.5 in 
the EMM state. Near the zone-center the signal is sharp 
and resolution limited. One can easily notice that the 
intensity with increasing energy shrinks and the width 
starts to increase. Following the sharp excitations with 
increasing energy, toward the zone boundary the signal is 
hampered, similar to observations in the optimal doped 
compound x=0.50. 



11 



The observation of nearly identical exchange parame- 
ters in Ndi-a^Sra^MnOa with x=0.49 and x=0.50 again 
emphasizes an identical FMM state in x=0.49 and 
x=0.50. The FMM state in Ndi_a,Sra,Mn03 seems to 
be less sensitive to the doping than the COSO phase 
which is severely damaged by the small amount of holes 
in Ndo.5iSro.49Mn03. 



VI. SUMMARY 

In conclusion, we have studied the microscopic origin of 
the COSO state in Ndo.sSro.sMnOs and the FMM state 
in Ndi-a^Sr^^MnOs with x=0.50 and x=0.49 by inelas- 
tic neutron scattering. Despite the complications arising 
from the twinning in these perovskite manganites with 
pseudo-cubic structure, we were able to observe excita- 
tions along the main-symmetry directions up to E?^12 
meV in the AFM ground state of x=0.5. We clearly find 
an anisotropic dispersion between the directions parallel 
and perpendicular to the orbital zigzag chains. Theo- 
retical descriptions of the spin-wave dispersion by the 
AFM CE-type ordering perfectly describe our experimen- 
tal data yielding results similar to the single-layered com- 
pound Lao.5Sri.5Mn04 (Ref. 15). The good description 
of the experimental dispersion gives strong support for 
the qualitative CE-type Goodenough model. A model 
with bond-centered charge order is even unable to explain 
the sign of the anisotropy. Although the Zener-polaron 
picture cannot explain the data, a superposition of both, 
Goodenough model and Zener-polaron as suggested in 



Ref. 14 can not fully be excluded. The stacking of mag- 
netic order along the c-direction, i.e. perpendicular to 
the planes of the zigzag chains, is rather strong. This 
finding can be understood by the fiattening of the crys- 
tal structure and by the fact this stacking arises from 
strong Mn^+-Mn^+ and Mn'^+-Mn'^+ interactions. 

The FMM state has been investigated in the half- 
doped compound Ndo.sSro.sMnOs and in the slightly 
electron-doped system Ndo.5iSro.49Mn03. Both concen- 
trations exhibit nearly the same spin- wave spectra, prov- 
ing that the FMM state is little affected by the slight 
difference in doping. In Ndo.sSro.sMnOs the magnon- 
dispersion along the main symmetry-directions [1,0,0], 
[0,1,1] and [1,1,1] can be observed. The FMM state 
seems to consist of two parts; (1) a spin- wave dispersion 
which can be easily described with a Heisenberg model 
taking only nearest-neighbors interactions into account 
and (2) weak additional intensities near the zone bound- 
aries over a wide energy band. The spin-wave dispersion 
in Ndo.5iSro.49Mn03 along the [1,0,0] and [0,1,1] direc- 
tions exhibits identical results with a steep dispersion and 
damped excitations near the zone-boundary. 
AcknoAvledgements 



This work was supported by the Deutsche Forschungs- 
gemeinschaft through the Sonderforschungsbereich 608. 
We thank S. Heijligen for the SQUID measurements. 
F. Kriiger acknowledges financial support from EPSRC 
under grant code EP/I004831/1. A. Nugroho thanks 
Cologne University and SFB 608 for hospitality. 



* Electronic address: ulbricli@ph2.uni-koelii.de 

^ Electronic address: frank . krugerQst-andrews .ac.uk 

■'■ Electronic address: braden@ph2.uni-koeln.de 

^ G. H. Jonker and J. H. V. Santen, Physica 16, 337 (1950). 

'^ J. H. V. Santen and G. H. Jonker, Physica 16, 599 (1950). 

^ Y. Tokura and Y. Tomioka, J. Magn. Mater. 200, 1 (1999). 

^ P. G. Radaelli, D. E. Cox, M. Marezio, and S.-W. Cheong, 

Phys. Rev. B 55, 3015 (1997). 
^ H. Kawano, R. Kajimoto, H. Yoshizawa, J. Fernandez- 

Baca, Y. Tomioka, H. Kuwahara, and Y. Tokura, Physica 

289 (1998). 
^ J. van den Brink, G. Khaliullin, and D. Khomskii, Phys. 

Rev. Lett. 83, 5118 (1999). 
^ S. Yunoki, T. Hotta, and E. Dagotto, Phys. Rev. Lett. 84, 

3714 (2000). 
^ R. Kajimoto, H. Yoshizawa, Y. Tomioka, and Y. Tokura, 

Phys. Rev. B 66, 180402 (2002). 
^ E. O. Wollan and W. C. Koehler, Phys. Rev. 100, 545 

(1955). 
^° J. Goodenough, Phys. Rev. 200, 564 (1955). 
^^ Z. Jirak, F. Damay, M. Hervieu, C. Martin, B. Raveau, 

G. Andre, and F. Bouree, Phys. Rev. B 61, 1181 (2000). 
^^ J. Herrero- Martin, J. Garcia, G. Subias, J. Blasco, and 

M. C. Sanchez, Phys. Rev. B 70, 024408 (2004). 
^^ A. Doud-Aladine, J. Rodriguez-Carvajal, L.Pinsard- 



Gaudart, M. Fernandez- Diaz, and A. Revcolevschi, Phys. 

Rev. Lett. 89, 097205 (2002). 
^^ D. Efremov, J. van den Brink, and D. Khomskii, Nature 

Mater. 3, 853 (2004). 
^^ D. Senff, F. Kriiger, S. Scheidl, M. Benomar, Y. Sidis, 

F. Demmel, and M. Braden, Phys. Rev. Lett. 96, 257201 

(2006). 
^^ G. Jung, V. Markovich, C. J. van der Beek, D. Mogilyan- 

sky, and Y. M. Mukovskii, Phys. Rev. B 72, 134412 (2005). 
^^ P. Rentier, O. Friedt, B. Bchner, M. Braden, and 

A. Revcolevschi, Journal of Crystal Growth 249, 222 

(2003). 
^^ D. Senff, P. Rentier, M. Braden, O. Friedt, D. Bruns, 

A. Cousson, F. Bouree, M. Merz, B. Biichner, and 

A. Revcolevschi, Phys. Rev. B 71, 024425 (2005). 
^^ D. Senff, O. Schumann, M. Benomar, M. Kriener, 

T. Lorenz, Y. Sidis, K. Habicht, P. Link, and M. Braden, 

Phys. Rev. B 77, 184413 (2008). 
2° H. Ulbrich, D. Senff, P. Steffens, O. J. Schumann, Y. Sidis, 

P. Rentier, A. Revcolevschi, and M. Braden, Phys. Rev. 

Lett. 106, 157201 (2011). 
^^ Y. Tokura and N. Nagaosa, Science 288, 462 (2000). 
^^ M. Tokunaga, N. Miura, Y. Moritoma, and Y. Tokura, 

Phys. Rev. B 59, 11151 (1999). 
^^ H. Yoshizawa, H. Kawano, J. Fernandez-Baca, H. Kuwa- 



12 



hara, and Y. Tokura, Phys. Rev. B 58 (1998). 
^^ J. A. Fernandez-Baca, P. Dai, H. Y. Hwang, C. Kloc, and 

S.-W. Cheong, Phys. Rev. Lett. 80, 4012 (1998). 
^^ H. Kawano-Furukawa, R. Kajimoto, H. Yoshizawa, 

Y. Tomioka, H. Kuwahara, and Y. Tokura, Phys. Rev. B 

67, 174422 (2003). 
^^ Y. Endoh, H. Hiraka, Y. Tomioka, Y. Tokura, N. Nagaosa, 

and T. Fujiwara, Phys. Rev. Lett. 94, 017206 (2005). 
^^ M. Hennion, F. Moussa, P. Lehouelleur, F. Wang, 

A. Ivanov, Y. M. Mukovskii, and D. Shulyatev, Phys. Rev. 

Lett. 94, 057006 (2005). 
^^ B. Kirby, J. Rhyne, H. Kaiser, H. Kuwahara, and 

Y. Tokura, J. Magn. Magn. Mater. 302, 237 (2006). 
^^ F. Ye, P. Dai, J. A. Fernandez-Baca, H. Sha, J. W. 

Lynn, H. Kawano-Furukawa, Y. Tomioka, Y. Tokura, and 

J. Zhang, Phys. Rev. Lett. 96, 047204 (2006). 
^° F. Ye, P. Dai, J. A. Fernandez-Baca, D. T. Adroja, T. G. 

Perring, Y. Tomioka, and Y. Tokura, Phys. Rev. B 75, 

144408 (2007). 
^^ S. Petit, M. Hennion, F. Moussa, D. Lamago, A. Ivanov, 

Y. M. Mukovskii, and D. Shulyatev, Phys. Rev. Lett. 102, 

207201 (2009). 
^^ R. Kajimoto, H. Yoshizawa, H. Kawano, H. Kuwahara, 

Y. Tokura, K. Ohoyama, and M. Ohashi, Phys. Rev. B 

60, 9506 (1999). 
^^ Y. Tokura, Rep. Prog. Phys. 69, 797 (2006). 
^^ H. Kuwahara, Y. Tomioka, A. Asamitsu, Y. Moritomo, 

and Y. Tokura, Science 270, 961 (1995). 
^^ H. Kawano, R. Kajimoto, H. Yoshizawa, Y. Tomioka, 

H. Kuwahara, and Y. Tokura, Phys. Rev. Lett. 78, 4253 

(1997). 
^^ J. Rodriguez-Carvajal, Physica B 192, 55 (1993). 
^^ In space group Imma it would align along the orthorhombic 

b-direction. 
^^ D. Pet it grand, H. Casalta, P. Bourges, and A. Ivanov, J., 

Physica B 806, 234 (1997). 



^^ V. V. Krishnamurthy, J. L. Robertson, R. S. Fishman, 

M. D. Lumsden, and J. F. Mitchell, Phys. Rev. B 73, 

060404 (2006). 
^° L. F. Feiner and A. M. Oles, Phys. Rev. B 59, 3295 (1999). 
^^ T. Hotta, S. Yunoki, M. Mayr, and E. Dagotto, Phys. Rev. 

B 60, R15009 (1999). 
^2 B. J. Sternlieb, J. P. Hill, U. C. Wildgruber, G. M. Luke, 

B. Nachumi, Y. Moritomo, and Y. Tokura, Phys. Rev. 

Lett. 76, 2169 (1996). 
^^ N. Shannon and A. V. Chubukov, Phys. Rev. B 65, 104418 

(2002). 
^^ W. Lv, F. Kriiger, and P. Phillips, Phys. Rev. B 82, 045125 

(2010). 
^^ C. I. Ventura and B. Alascio, Phys. Rev. B 68, R020404 

(2003). 
^^ O. Sikora and A. M. Oles, Physica B 359-361, 1300 (2005). 
^'^ L. E. Gontchar and A. E. Nikiforov, Phys. Rev. B 66, 

014437 (2002). 
^^ L. Jongh, Magnetic properties of layered transition metal 

compounds (Kluwer Academic, 1990). 
^9 T. G. Perring, G. Aepph, S. M. Hayden, S. A. Carter, 

J. P. Remeika, and S.-W. Cheong, Phys. Rev. Lett. 77, 

711 (1996). 
^° H. Hwang, P. Dai, S. Cheong, G. Aepph, D. Tennant, and 

H. Mook, Phys. Rev. Lett. 80, 1316 (1998). 
^^ J. Zhang, F. Ye, H. Sha, P. Dai, J. A. Fernandez-Baca, and 

E. W. Plummer, Journal of Physics: Condensed Matter 

19, 315204 (2007). 
^^ A. Zheludev, ResLib 3.3 (Oak Ridge National Laboratory, 

Oak Ridge, 2006, 2006). 
^^ M. C. Martin, G. Shirane, Y. Endoh, K. Hirota, Y. Morit- 
omo, and Y. Tokura, Phys. Rev. B 53, 14285 (1996). 
^^ In order to compare the J's in the FMM state with those 

in the CE-type phase one needs to divide the latter by a 

factor two. 



